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ABSTRACT 

A  FORTRAN  IV  computer  program  was  employed  to  conduct  a 
statistical  analysis  of  data  collected  during  fleet  anti- 
submarine warfare  exercises.   The  object  of  the  investigation 
was  the  identification  of  those  variables  which  had  greatest 
influence   on  a  destroyer's  ability  to  detect  a  submarine 
under  certain  conditions. 

The  variables  were  treated  as  a  random  vector  arising 
from  one  of  two  multivariate  normal  populations  with  common 
covariance  matrix.   An  artificial  regression  relation  was 
formulated  to  facilitate  development  of  a  linear  discriminant 
function  in  a  subset  of  those  variables  found  to  be  of  dominant 
importance.   This  latter  subset  was  identified  by  examination 
of  multiple  correlation  coefficients. 

The  discriminant  function  was  found  to  be  seventy  five 
per  cent  effective  in  classifying  the  experimental  data 
correctly. 
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I.   INTRODUCTION 

The  sonar-equipped  destroyer  is  the  primary  surface 
vessel  employed  by  the  U.S.  Navy  in  antisubmarine  warfare. 
Detection  of  the  submarine  necessarily  precedes  any  offen- 
sive action  the  destroyer  can  take  against  it.   It  was  the 
objective  of  the  author  to  identify  those  variables  which 
had  the  greatest  influence  on  the  effectiveness  of  a 
destroyer  in  detecting  a  submerged  submarine  during  fleet 
exercises.   Consideration  was  limited  to  detection  capabil- 
ities while  using  sonar  equipment  in  the  active  mode. 

Reference  2  lists  many  of  the  large  number  of  variables 
which  influence  the  effectiveness  of  a  destroyer  in  submarine 
detection.   Certain  of  these  variables  can  be  regarded  as 
control  variables  in  the  sense  that  the  value  of  those  vari- 
ables can  be  chosen  and  fixed  by  the  destroyer  commander  with 
a  certain  degree  of  freedom.   For  example,  the  destroyer's 
speed  may  range  from  zero  to  thirty  knots.   Furthermore,  a 
particular  destroyer  can  be  equipped  with  one  of  several  types 
of  sonar  system,  and  can  be  manned  with  differing  numbers  of 
sonar  technicians  of  various  skill  levels.   These  latter 
variables  serve  as  examples  of  control  variables  whose  value 
can  be  adjusted  by  agencies  exogenous  to  the  ship.   Similar 
comments  can  be  made  concerning  the  submarine  involved  in  an 
encounter  with  a  destroyer.   No  attempt  was  made  by  the  author 


to  consider  any  variable  as  a  control  variable;  all  were 
regarded  as  completely  random  variables. 

Based  on  personal  experience  in  the  Destroyer  Force,  the 
author  chose  a  set  of  twenty-two  variables  whose  examination 
was  expected  to  be  most  fruitful  for  the  purpose  stated. 
Four  groups  of  variables  were  considered:   environmental 
factors,  destroyer  characteristics,  submarine  characteristics, 
and  tactical  factors.  "Table  1  lists  the  variables  included 
in  each  group.   Qualitative  variables,  such  as  destroyer  hull 
type  were  assigned  numerical  codes.   All  values  were  recorded 
at  the  time  detection  was  achieved,  or  -  in  those  cases  where 
detection  did  not  occur  -  when  the  submarine  had  approached 
to  a  specified  range.   This  range  was  specified  by  the  U.S. 
Navy  Antisubmarine  Warfare  Data  Center  (NADAC)  as  a  standard 
for  each  sonar  system.   The  standards  are  listed  in  Ref.  2. 


Table  I.   Variables  by  Group 


Environmental 

Destroyer 

Submarine 

Tactical 

Geographic 

Hull  type 

Hull  type 

Relative  speed 

location 

between 
opponents 

Present  weather 

Hull  number 

Hull  number 

Distance  between 
opponents 

Sonic  layer 

Sonar  type 

Depth 

Aspect  angle  of 

depth 

submarine 

Wind  speed 

Speed 

Speed 

Relative  bearing 
from  destroyer 
to  submarine 

Sea  state 

Self  noise 

Surface  turbu- 

Technician com- 

lence 

plement 

Bathythermo- 

graph profile 

type 

Water  depth 

The  data  examined  was  supplied  by  NADAC  on  request,  having 
been  collected  during  sixteen  fleet  exercises  conducted  in 
the  Atlantic  and  Pacific  Oceans  during  the  period  1962  to 
1966,  inclusive.   The  destroyer  types  considered  included  DLG, 
DL,  DDG,  DDR,  DD ,  DEG  and  DE .   Each  destroyer  was  equipped 
with  one  of  the  following  AN/SQS  sonar  systems:   23,  26,  29, 
or  30.   Submarine  types  encountered  included  all  types  in  the 
U.S.  fleet  during  this  period,  ranging  from  the  World  War  II 
vintage  fleet  submarine  to  the  nuclear-powered  ballistic 
missile  submarine. 

Each  observation  of  the  twenty-two  variables  was  regarded 
as  a  random  vector,  x,  arising  from  one  of  two  multivariate 
normal  populations  with  common  covariance  matrix:  N (m, ,  V) , 
or  NCnu,  V).   The  former  distribution  representing  an  encounter 
which  resulted  in  detection  of  the  submarine  by  the  destroyer, 
and  the  latter  distribution  applying  otherwise. 

Maximum  likelihood  estimators  were  used  for  the  mean 
vectors,  covariance  matrix,  and  all  other  parameters  employed 
in  the  analysis. 

a  FORTRAN  IV  computer  program  was  prepared  for  analysis 
of  the  data  and  was  utilized  with  the  U.S.  Naval  Postgraduate 
School  IBM-360  computer.   A  discussion  of  the  program,  a  flow 
chart,  and  the  program  appear  below. 

In  all,  observations  from  ninety  two  encounters  were 
included  in  the  analysis.   Twenty-seven  encounters  resulted  in 
submarine  detection  by  the  destroyer,  and  sixty-five  did  not. 


II.   DISCUSSION 

A.   PRINCIPAL  STEPS  OF  THE  ANALYSIS 

The  following  discussion  is  concarned  with  the  principal 
concepts  employed  in  the  analysis ,  without  regard  to  the 
computer  programming  techniques  employed.   A  full  discussion 
of  the  program  appears  later. 

The  linear  discriminant  function 

x1  .   V_1  (n^  -  m2)  "  ^(^  +  m2)  '  V-1  (m,  -  m_2)  =  0    (1) 

is  related  to  the  artificial  regression  proposed  by  Fisher 
in  Ref.  5 

y.    =  b'  (x_.  -  c)  .  (2) 

In  the  foregoing,  x.  is  the  observation  vector  with  j  =  1, 

2,  ...  ,  n,  +  n~;  there  are  n,  observations  in  the  first 

population  with  parameters  (m, ,  V)  and  n~  observations  in 

the  second  population.   The  artificial  dependent  variable, 

y.,  takes  the  value  n_/(nn  +  n~)  if  x.  is  from  population 
Jj  2    12      —j  r    c 

one,  and  the  value  -n./Cn,  +  n„)  otherwise.  The  vector  c 
is  a  constant  vector  with  each  element  being  the  weighted 
grand  mean  of  observations,  i.e.,  (n,  x.  +  n„  x?) / (n,  +  n2) . 

Equation  (2)  yields  artificial  regression  coefficients, 
b,  which  are  proportional  to  the  discriminant  function  coef- 
ficients, V~   (m,  -  m  ) ,  of  equation  (1).   (This  may  be 
verified  by  solution  of  the  normal  equations  in  the  components 
of  b.) 


Based  on  the  above  relationship  between  the  discriminant 
and  regression  functions,  equation  (2)  was  utilized  to 
determine  a  subset  of  variables  which  could  be  regarded  as 
being  of  principal  importance  in  the  discriminant  function. 
The  technique  employed  is  discussed  in  Ref.  7.   A  multiple 
correlation  coefficient,  R,  ,  ,  was  calculated  from  equation 
(1)  using  the  full  set  of  twenty-two  independent  variables. 
Then,  each  independent  variable,  x.,  was  eliminated  in 
iterative  fashion  and  a  corresponding  multiple  correlation 
coefficient,  R~  . ,  was  calculated  based  on  the  remaining  twenty- 
one  independent  variables.   The  significance  of  each  of  the 
twenty-two  variables  in  the  regression  was  then  tested 
according  to  the  following  scheme.   If  the  following  relation- 
ship was  satisfied,  then  the  deleted  variable,  x.,  was  discard- 
ed as  contributing  little  to  the  regression: 

q(Rll  "  R2i)/(1  "  Rll} (pl  "  P2}  -  F   *  (2a) 

In  this  relationship  p1  was  the  number  of  variables  in  the 
original  regression  (22) ,  p_  was  the  number  in  the  reduced 
set  (21)  ,  F  was  the  F  oc-  variate  with  (p,  -  p„)  and  q  degrees 
of  freedom,  respectively;  q  had  the  value  (n,  +  n„)  - 

(p1   +    1) • 

The  above  technique  revealed  that  all  but  six  of  the 
original  variables  contributed  little  to  the  regression,  and  - 
therefore  -  to  the  discriminant  function.   The  value  of  R,, 
was  found  to  be  0.646  for  the  original  variable  set  and  had 
the  value  0.528  for  the  six  significant  variables.   The  latter 


value  was  calculated  from  a  regression  based  on  these  six 
variables  alone.   That  is,  with  a  seventy-three  per  cent 
reduction  in  the  number  of  independent  variables,  the  corre- 
lation coefficient  was  reduced  by  only  eighteen  per  cent. 
Table  II  lists  the  six  significant  variables  in  order  of 
importance,  the  value  of  the  correlation  coefficient  for  the 
set  of  twenty-one  variables  remaining  upon  deletion  of  each 
of  the  significant  variables,  and  the  reduction  in  the  multiple 
correlation  coefficient  which  resulted  upon  deletion  of  the 
variable . 

With  the  six  significant  variables  thus  identified,  a 
linear  discriminant  function  was  constructed,  and  its  properties 
tested.   The  specific  form  of  equation  (1)  used  at  this  point 
was 

x'.  d  -  -k  =  0,  (3) 

—j    —    2 

4-Vi 

where  x.  is  the  j    6X1  observation  vector,  d  is  the  vector 
-D  - 

of  discriminant  coefficients  and  k  is  a  constant.   The  relation 
between  equations  (1)  and  (3)  is  given  by 

d  =  V    (m,  -  m„) 

k  =  (m,  +  m  )  *  V    (m,  -  mj 

With   the  value  of  all  parameters  calculated,  equation  (3) 
took  the  following  form: 

-0.165x,  +  0.015xo  +  0.411xo  -  0.020x.  +  0.506X,. 
1  2  3  4  o 


(4) 


-  0.008x,  >  8.28  . 
o  — 


Table  II.   Significant  Variables. 


_      Per  cent  Reduction 
i 2i in  Correlation 

5  Type  of  weather  0.565  12.5 

2  Submarine  hull  number     0.604  6.5 

6  Sonic  layer  depth         0.60  6  6.2 
1      Geographic  location       0.60  8             5.9 

3  Submarine  hull  type       0.612  5.3 

4  Submarine  depth  0.617  4.5 

An  observation  would  be  classified  as  arising  from  population 
one  if  the  relation  was  satisfied,  and  as  arising  from  popula- 
tion two  otherwise.   With  the  formulation  of  equation  (4),  the 
observations  were  then  to  be  classified  using  the  equation, 
and  the  results  achieved  to  be  considered  a  measure  of  the 
validity  of  the  discriminant  function. 

Anderson  discusses  in  Ref.  1  the  probability  of  mis- 
classification  associated  with  a  linear  discriminant  function. 
Under  the  assumption  that  the  cost  of  an  error  of  misclass- 
ification  is  equal  for  each  of  the  two  possible  errors,  the 
probability  of  misclassif ication,  P (M) ,  is  given  by 


P  (M)  =  P  (Y  >  /a/2'   ) 

where  Y  is  a  standard  normal  variate  and  a  takes  the  value  of 
(m,  -  nu)  V    (m,  -  m„) .   From  the  above  equation,  it  was 
found  that  the  value  of  P (M)  associated  with  equation  (4)  is 
0.18. 

Upon  classifying  the  ninety-two  observations  of  the  six 
significant  variables  -  using  equation  (4)  -  it  was  found  that 
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eight  (30%)  of  the  twenty-seven  cases  where  detection  occurred 
were  misclassif ied  as  cases  of  non-detection.   Similarly, 
fourteen  (22%)  of  the  sixty-five  cases  of  non-detection  were 
misclassif ied.   Thus,  the  overall  error  of  misclassif ication 
was  twenty-five  per  cent. 

In  view  of  the  small  sample  sizes  available,  particularly 
in  the  case  of  population  one  (twenty-seven  observations) ,  the 
difference  between  the  theoretical  error  of  misclassification 
and  the  observed  value  was  considered  acceptable.   The  fact 
that  the  error  of  misclassification  in  each  population 
decreased  with  sample  size  was  viewed  as  an  encouraging  sign. 

Consequsntly ,  it  was  concluded  that  the  six  variables 
cited  above  were,  in  fact,  of  dominant  importance  in  destroyer  - 
submarine  encounters . 

It  was  noted  with  interest  that  the  character  of  these 
variables  is  such  that  the  destroyer  commander  has  virtually 
no  means  available  to  influence  them  to  his  advantage.   There 
seems  to  be  little  direct  action  that  could  have  been  taken 
by  the  destroyer  commander  to  increase  significantly  the  like- 
lihood of  detecting  his  submerged  opponent. 

On  the  other  hand,  the  submarine  commander  enjoys  a  distinct 
advantage.   If  his  depth  and  location  are  chosen  judiciously, 
and  his  submarine  is  of  the  preferred  design,  he  may  proceed 
with  confidence  of  not  being  detected  by  his  surface  adversary. 

B.   THE  COMPUTER  PROGRAM 

The  computer  program  employed  in  the  analysis  contained 
four  principal  steps.   First,  the  data  were  read  into  storage, 
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constants  were  defined,  and  the  data  were  transformed  into  a 
form  to  facilitate  subsequent  calculations.   Second,  the 
multiple  correlation  coefficient  for  the  twenty-two  independent 
variables  was  calculated.   Third,  using  the  scheme  described 
above,  the  significant  variables  were  identified.   Finally,  the 
discriminant  function  was  formulated  and  tested.   The  discus- 
sion below  is  concerned  with  the  mechanics  and  logic  of  the 
entire  program.   A  flow  chart  of  the  program  is  included  as 
Appendix  A,  and  is  followed  by  a  listing  of  the  program  itself. 

The  following  subprograms  which  were  prepared  by  the 
computer  manufacturer  formed  a  substantial  part  of  the  program: 
CORRE,  ORDER,  MINV,  MULTR,  DMATX .   A  full  discussion  of  these 
subprograms  may  be  found  in  Ref.  6.  They  are  discussed  below 
briefly  where  necessary  to  maintain  continuity. 

After  storage  requirements  and  constants  are  defined  in  the 
main  program,  control  is  passed  to  the  subroutine  DATA.   Here, 
a  DC  loop  is  used  to  read  the  raw  data  into  storage  as  a 
three-dimensional  array  X(I,J,K).   (The  indices  of  X  are 
observation  number,  variable  number,  and  population  number 
respectively.)   As  each  card  of  data  is  read,  the  values  of 
four  variables  are  calculated  and  the  data  are  transformed  to 
a  one  dimensional  array  XD.   The  latter  step  is  required  by 
the  subsequent  use  of  the  manufacturer-supplied  subprograms 
mention  above . 

The  first  value  calculated  is  that  of  the  dependent  var- 
iable y.  which  is  stored  in  X(I,23,K).   Next  the  relative 
bearing  from  destroyer  to  submarine  is  converted   to 
radian  measure  using  the  simple  subprogram  REL.   With  destroyer 
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course  and  submarine  true  bearing  as  input  variables,  REL 
calculates  the  relative  bearing  as  the  output  variable  RB. 
This  value  is  then  converted  to  radian  measure  and  stored  in 
X(I,12rK).   Then,  relative  speed  is  computed  from  destroyer 
speed,  submarine  speed,  and  the  angle  between  the  two  vessels' 
courses.   This  value  is  stored  in  X(I,5,K).   The  value  of 
sonic  layer  depth  is  stored  in  X(I,22,K)  based  on  the  fleet 
exercise  from  which  the  observation  was  drawn.   Finally,  the 
destroyer's  manning  level  in  the  sonar  technician  (ST)  rating 
is  used  to  calculate  a  parameter  reflecting  the  rate  distri- 
bution.  This  is  accomplished  by  calculating  a  sum  of  weighted 
values  and  storing  the  result  in  X(I,13,K).   For  each  destroyer, 
the  number  of  chief  petty  officers  is  multiplied  by  six,  first 
class  petty  officers  by  five,  and  so  on,  until  a  product  is 
formed  for  all  skill  levels  in  the  ST  rating.   The  sum  of  these 
products  is  then  taken  as  the  measure  of  manning  level  for  that 
destroyer. 

With  the  data  transformation  completed,  control  returns 
to  the  main  program  and  the  multiple  linear  regression 
equations  are  formulated.   The  first  subprogram  called  is 
CORRE.   Given  as  input  variables  the  number  of  observations 
NS,  the  number  of  variables  M,  and  the  one  dimensional  array 
of  observations,  XD  -  CORRE  computes  the  vector  of  means 
CMEAN,  the  vector  of  standard  deviations  PS,  the  MxM  matrix  of 
sums  of  cross-products  of  deviations  from  means  D,  and  the 
upper  triangular  portion  of  the  MxM  symmetric  matrix  of 
correlation  coefficients  R. 
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After  the  subscripts  of  the  variables  are  stored  in  the 
array  ISAVE  in  the  order  in  which  the  variables  are  arranged 
in  the  array  XD,  control  is  passed  to  subroutine  ORDER. 

Given  the  number  of  variables,  the  matrix  of  correlation 
coefficients/  the  subscript  number  of  the  dependent  variable, 
and  the  array  ISAVE,  subroutine  ORDER  computes  the  22x22  matrix 
of  correlation  coefficients  among  independent  variables  R,  , 
and  the  vector  of  correlation  coefficients  between  the  dependent 
variable  and  the  independent  variables  R, „ . 

Next,  the  subroutine  MINV  is  called  to  calculate  the 
inverse  of  R, , .   This  matrix  plus  the  vector  of  means,  the 
vector  of  standard  deviations  and  the  diagonal  elements  of 
the  matrix  D  are  then  passed  to  the  subroutine  MULTR.   This 
subprogram  computes  the  vector  of  regression  coefficients  W 
and  the  multiple  correlation  coefficient  ANS(2).   (The  vector 
W  and  ANS(2)  correspond  to  the  vector  b  of  equation  (2)  and 
the  variable  R, ,  of  the  relation  (2a)  discussed  earlier.) 

A  DO  loop  is  then  used  to  eliminate  each  independent 
variable  from  the  regression,  and  twenty-two  multiple  corela- 
tion  coefficients  for  each  combination  of  twenty-one  variables 
is  computed.   The  test  statistic  of  relation  (2a)  is  computed 
for  each  variable  x.  and  stored  in  C(J) .   As  each  new 
multiple  correlation  coefficient  is  calculated,  the  value  of 
the  test  statistic  is  stored  temporarily  in  TEM3  for  each 
test.   This  value  is  then  compared  to  the  value  of  the  F(l,69) 
variate  at  the  5%  significance  level.   If  the  variable  with 
subscript  J  is  found  to  be  of  little  significance,  then  a  zero 
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is  stored  in  array  MK,  otherwise  MK(J)  will  contain  the  value 
of  the  subscript. 

On  completion  of  the  above,  the  values  of  the  test  sta- 
tistic are  arranged  in  order  of  magnitude  in  the  array  C,  and 
deleted  from  the  array  I SAVE. 

At  this  point  the  number  of  significant  variables  is 
checked.   This  step  is  necessitated  by  the  fact  that  the  array 
TST  contains  values  of  the  F(f, ,f-)  statistic  with  f1  =  1  and 
6  _<  f~  _<  14.   The  statistic  is  used  later  in  the  program  and 
the  results  of  the  program  would  be  invalid  if  the  number  of 
significant  variables  were  not  within  the  range  given  for  f„. 
(The  range  chosen  was  based  on  preliminary  trials  with  the 
program.   It  was  expected  that  there  would  be  six  significant 
variables . ) 

The  subprograms  ORDER,  MINV  and  MULTR  are  employed  as  before 
with  the  regression  equations  being  based  on  the  significant 
variables  only.   Deletion  of  the  indices  of  the  insignificant 
variables  from  the  array  ISAVE  results  in  the  appropriate 
adjustments  in  the  subprograms.   The  test  statistic  TEM3  is 
compared  to  the  appropriate  F  variate  TST (IT),  and  the  index 
is  incremented.   If  the  test  is  not  satisfied,  i.e.,  if  the 
discarded  variables  have  significance  as  a  group  (an  unexpected 
result) ,  a  value  of  one  is  stored  in  ICK  and  this  will  result 
in  a  change  of  the  flow  of  the  program  later. 

A  DO  loop  follows  the  above  and  rearranges  the  order  of 
the  variables  in  the  array  XD  in  preparation  for  calculation 
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of  the  discriminant  function.   Control  then  passes  to  the 
subroutine  DMATX. 

Given  the  number  of  (significant)  variables  L,  the  number 
of  populations  LX,  the  number  of  observations  in  each  popula- 
tion N,  and  the  array  of  data  XD  -  DMATX  calculates  the  vector 
of  means  XBAR  and  the  common  covariance  matrix  D.   The  inverse 
of  D  is  then  calculated  by  subroutine  MINV  and  control  passes 
to  subroutine  DISC.   This  latter  subroutine  is  a  simple 
FORTRAN  IV  equivalent  of  equation  (1) . 

Next,  each  observation  is  classified  using  the  discriminant 
function  just  formulated.   The  value  of  the  function  is  tested 
against  zero  and  the  observation  classified  accordingly.   For 
each  population,  the  value  of  KOUNT  is  incremented  each  time 
an  observation  is  improperly  classified. 

After  all  observations  have  been  classified,  the  value  of 
ICK  is  checked.   As  discussed  above,  if  ICK  has  value  one,  the 
previously  discarded  variables  have  group  significance.   If 
this  is  the  case,  a  warning  to  that  effect  is  printed,  the 
index  of  the  most  significant  deleted  variable  is  added  to  the 
array  MK,  control  returns  to  statement  number  45,  the  new 
group  of  significant  variables  are  tested  and  a  new  discriminant 
function  is  found.   If  ICK  has  a  value  of  zero,  the  discrim- 
inant function  coefficients  are  printed  and  the  program 
terminates . 
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III.   CONCLUSION 

Analysis  of  data  collected  during  fleet  exercises  conducted 
during  the  period  1962  to  1966  inclusive  strongly  suggests 
that  the  variables  listed  in  Table  II  had  a  dominant  influence 
on  destroyer  -  submarine  encounters.   Specifically,  the 
successful  detection  of  a  submarine  by  a  destroyer  utilizing 
a  sonar  of  the  type  specified  in  the  active  mode  was  influenced 
more  strongly  by  these  six  variables  than  by  the  remaining 
sixteen  considered. 

The  discriminant  function,  equation  (4) ,  was  75%  effective 
in  classifying  encounters  and  should  be  useful  as  a  predictive 
tool  in  assessing  the  likelihood  of  detection  under  known 
environmental  conditions  and  against  known  submarine  types. 

Consideration  of  the  variables  found  to  be  significant 
indicates  that  the  tactical  advantage  was  with  the  submarine 
commander  during  encounters  of  the  type  considered. 
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FLOW  CHART  OF  COMPUTER  PROGRAM 


READ   LX, 
M,    N 


JL 


ICK=0 


IT=0 


]f 


KI=M-1 


\L 


NS=N(1) 
+    N(s) 
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o 


DEFINE 
ARRAY 
TST 


CALL  DATA 


Y_ 


CALL  CORRE 


1  =  1 


V 


ISAVE (I) 
=  I 


I  =  1+1 


19 


0 


CALL  ORDER 


CALL  MINV 


CALL  MULTR 


!L 


WRITE  W 


/ 


jL 


WRITE 
LANS(2) 


TEM1 
ANS(2) **2 


20 


0 


© 


DN    = 
(l-TEMl)/69 

\ 

' 

KD=KI-1 

V 

' 

J=l 

11 

' 

L=0 

p 

_J 

/ 

JI=1 

21 


JI==JI+1 


L=L+1 

N 

' 

ISAVE(L) 

=    JI 

CALL  ORDER 


CALL  MINV 


22 


PRINT  MULT.  CORRELATION 

COEFFICIENT  FOR  JTH 

GROUP  OF  21  VARIABLES 


TEM2 
=ANS (2) **2 


C(J)=TEM3 


27 


-> 


MK(J)=J 


MK(J)=0 


30 


J=J+1 


< 


e 


KJ=1 


A    ) — — — > 


STR=C(KJ) 


IND=KJ 


v 

0 


24 


0 


KL=1 


KL=KL+1  U- 


IDX(KJ)=IND£ 


-I 


C(IND)=0 


KJ=KJ+1 


STR=C(KL) 


y_ 


IND=KL 


25 


I SAVE (L) 
=MK(JJ) 


JJ=JJ+1 


26 


WRITE  ISAVE  (INDICES  OF 
SIGNIFICANT  VARIABLES) 


CALL  MINV 


CALL  MULTR 


WRITE 
ANS  (2) 


27 


© 


TEM2  = 
ANS  (2) **2 


TEM3= (TEM1- 

TEM2)/ 
(DN*(KI-L) ) 


IT=IT+1 


LK=0 


K=l 


-> 


ICK=1 


28 


© 


© 


NK=N(K) 


JK=1 


LK=LK+1 


1. 


N2=LK-NK 


IE=1 


& 


JR   = 
I SAVE (IE) 


29 


© 


N2=N2+NK 


XD(N2)  = 
X(JK,JR,K) 


IE==IE+1 


LK=N2 


-> 


JK=JK+1 


K=K+1 


30 


© 


CALL  DMATX 


WRITE  XBAR(II) ,  XBAR(JX) 
(VECTORS  OF  MEANS) 


CALL  MINV 


1 


CALL  DISC 


KN=1 


v 


KOUNT=0 


Q 


31 


0 


N1=N(KN) 


IL=1 


SUM=0 


«:- 


<D 


32 


© 


IG  = 
ISAVE(JL) 


IND=IND+1 


X 


SUM=SUM+ 
(IND) *X(IL, 
IG,KN) 


JL=JL+1 


/SUM-\ 

NN= 

-2 

\C(L+1)/ 

*>> 

<o 

* 

1° 

NN=1 

^ 
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33 


0 


& 


WRITE 
KOUNT,KN 


WRITE 
WARNIN( 


f 


KOUNT= 

KOUNT+1 


IL=IL+1 


KN=KN+1 


LK=L+1 


-0 


I, 


-■© 


34 


© 


© 


IND=IDX(L+1 


MK(IND)=IND 


ICK=0 


± 


WRITE  C,  VECTOR  OF 
DISCRIMINANT  FUNCTION 
COEFFICIENTS 


35 


COMPUTER  PROGRAM 

DIMENSION    X(    65,23,2)  ,X3AR(44)  ,     D(529)t     CMEAN(23)f 
/XD(  2116 ) ,PS(23)  ,h(2 ) ,R( 523)  ,B<?3)  ,W(23)  ,T(23)  , 
/ISAVE(23), Fill <^64) , R12( 22 ), SW<2  3) ,ANS( 1C) , MK ( 22) , 
/IDX(22)  ,TST(10)  ,C(23) 

READ(5,1)     LX,M,N 

ICK  =  C 

IT=0 

KI=M-1 

NS=N(1 )+N(2) 

DATA    T ST /2. 36, 2. 39, 2. 42, 2. 45, 2. 49, 2. 5 5, 2. 60, 2. 66, 
/2.8C2.Q5/ 
C 

C    TRANSFORM  RAW  DATA  FOR  MULTIPLE  REGRESSION  ANALYSIS 
C 

CALL  DATA< M,N,NS,X,xn) 

CALL  COR  RE ( MS , M , 1 , XD,C MEAN , PS ,D,R,B,W,T) 

DC)  1C  1=1,  K I 
10  ISAVE(  I  )=1 

CALL  ORDER (  A , R , M , K  I  , I SA VE , R 1 1 , R  1  2 ) 

CALL  MI  MM  Pll , KI ,DT,W, T  ) 

CALL  MULTMNS,KI  ,CMtAN,  PS,B,R11  ,R12  , 1  S  AVE  ,  W,  SW  ,T  ,  ANS  ) 

WRITE (6, EC C ) 

WRITER  ,5C1)  (W(  I  A)  ,  I  A=1,KI  ) 

WRITE(6,81C) 

WRIT E< 6 ,5011  AMS(2) 

TEM1=( ANSt ?) )**2. 

DN=( 1.-TEM1 )/69. 

KD=KI-1 

DO  30  J=1,KI 

L  =  0 

DO  20  J  1=1 f KI 

IF(  JI.ECUJ)  GO  TO  20 

L  =  L  +  1 

ISAVF  ( L  )=J  I 
20  CCNTINUE 

CALL  ORDER  (M,R,,M,KD,  ISAVE,R11,R12) 

CALL  MINV( Pll ,KDtDT,W,T ) 

CALL  MULTR ( N S , KD, CME A N , PS , B , R 1 1 ,P 1 2 , I S AVE , W, S W , T , ANS ) 

WRITE(6,?20)  J 

WPITE(6,5C1)  ANSt 2  I 

TEM2=(ANS(21 )**2. 

TEM3=(TEM1-TEM2 )/DN 

C(J)=TEM3 

IF(TEM3.GT.3.98)  GO  TO  27 
25  MK(J)=0 
C 

C     MK(J)=0  IMPLIES  THAT  X(J)  DOES  NOT  CONTRIBUTE 
C   SIGNIFICANTLY  TO  THE  CORRELATION  COEFFICIENT. 
C 

GO  TO  3C 
27  MK(J)=J 
30  CONTINUE 

DO  4C  KJ=1 , 22 

STR=C(KJ) 

IND=KJ 

DO  35  KL=1,22 

IF(C(KL).LE.STR)  GO  TO  35 

STP=C(KL) 

IND=KL 
35  CONTINUE 

IDX(KJ) =IND 
40  C<  IND)  =  C. 

C     IDX  CONTAINS  THE  INDICES  OF  THE  • ■ I NDEPENDENT  •  ■ 


36 


VARIABLES  ORDERED  WITH  RESPECT 
THE  VARIABLE'S  CONTRIBUTION  TO 
COEFFICIENT 


TO  SIGNIFICANCE  OF 
THE  CORRELATION 


45 


60 


C 
C 

c 

C 


L  =  0 

DO 

IF( 

L=L 

ISA 

CON 

IF( 

WRI 

WRI 

CAL 

CAL 

CAL 

WRI 

WRI 

TFM 

TEN 

IT  = 

IF( 


60  J 
MM  J 
+  1 

VE(  L 
TINU 
(L. 
TE(  6 
TEU 
L  OR 
L  MI 
L  MU 
TEU 
TE(6 
2=(A 
3=(T 
17  +  1 
TEM3 


J=1,KI 

JUEQ.O)  GO  TO  60 

)=MK( JJ) 

E 

LT.( ).0R. (L.GT. 15)  )  GO  TO  1010 

,830) 

,2)  (  ISAVE(  IV) ,  IV=1,L) 

DER(MfR,M,Li ISAVE, Rl 1 ,R12) 

NV( F 11,L,DT,W, T) 

LTR(NS,L,CMEAN,PS,B,Rll,R12,ISAVE,W,SW,TfANS) 

,840) 

,501 >  ANS< 2) 

NSC  2) >**2. 

EM1-TEM2 )/<DN" (KI-L) ) 


.GT.TST(IT)  )  ICK=1 


C 
C 
C 

C 


DISCRIMINANT  FUNCTION  IS  FORMED  BASED  ON 
SIGNIFICANT  «• INDEPENDENT • •  VARIABLES. 

LK  =  C 

DO  EO  K=l»2 

NK=N(K) 

DO  70  JK=1,NK 

LK^LK+1 

N2=LK-NK 

DC  7  0  IE=1 

JR=ISAVE(I 

N2=N2+NK 

XL(N2)=X( JK,JR,K) 

LK=N2 

CALL  DMATX( LX,L,N,XD,XBAR,D,T) 

WRITE (6,99 5) 

DC)  77  I  1=1  ,  L 

JX=I I+L 

WRIT E( 6, 999)  XBAR(I 

CALL  MINV( D,L, DT,T, 

CALL  01  SCt  XriAP.  ,L,D, 


THE 


,L 
E) 


70 
80 


77 


I), 

C) 
T,C 


XBAR( JX) 


INPUT  TATA 
DISCRIMINANT 


IS  CLASSIFIED 

FUNCTION 


ACCORDING    TO    THE    CURRENT 


79 


5 

88 
85 
81 

90 


1000 


x< 

4, 


DO    81-    KN  =  1  ,2 

KOUNT=0 

N1=N(KN) 

DO    85    IL=1 ,N1 

SUM=0. 

IND=C 

DO    79    JL=1 ,L 

IG=ISAVE(JL) 

IND=INP+1 

SUM=SUM+C(  IND) 

IF(SUM-CU. +1) ) 

NN=2 

GC    TO    F8 

NN=1 

IFINN.NE.KN) 

CONTINUE 

WRITE  (6,6C1)     KOUiMT, 

IF(  ICK. EC. C)     GO    TO 

WPITE(6,91 ) 

IND=IDX (L+l ) 

MK( I  NO )= IND 

ICK=0 

GO    TO    45 

LK=L+1 


IL,IG, 
5,5 


KN) 


K0UNT=K0UNT+1 


KN 
1C00 
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1010 
1015 

1 
2 

91 

500' 

501 

6C1 

800 

310 

820 

830' 

340 

995 

999 


WRITE(6 
WPITE(6 
STOP 
WRITHo 
FORMAT ( 
STOP 
FORMAT! 
FORMAT  ( 
FORMAT { 

'•SIGNIF 
FORMAT ( 
FORMAT! 
FORMAT ( 
FORMAT! 
FORMAT ( 
FORMAT  ( 

"GROUP  ' 
FORMAT! 

'•CONTR  I 
FORMAT ( 

"  VARIA8 

FORMAT ( 

FORMAT! 

EN 


,500) 
,501) 


( C( JS) , JS=1,LK) 


,  1 C  1  5  )     L 

•WARMING* , • L=» , 110) 


M4 
10X 
1CX 
ICA 
25X 
8E1 
1HC 
IPC 
1HC 
1HC 
,14 
1H0 
BUT 

1  HO- 
LES 

2  ex 

2E4 
D 


) 

,2 

,' 

NC 
,' 

5. 
,1 
,1 


RDED   VARIABLES  HAVE  GROUP  ', 
IMINANT  FUNCTION  COEFF IC  I  ENTS •  ) 


) 
,  • 

ir 
•  ) 
o. 


014) 

01  SCA 

E«  ) 

DISCR 

3) 

35,  4X, •MISCLASSIFICATIHNS  IN  GROUP* 

OX,'  S 

M  UL  T  I 

MULT  I 


TANOARD  REGRESSION  COEFFICIENTS' 

PLE  CORRELATION  COEFFICIENT') 

PLE  CORRELATION  COEFFICIENT  FOR  • 


16) 
) 


INDICES  OF  VARIABLES  MAKING  SIGNIFICANT  • 
N  TO  THE  CORRELATIJN  COEFFICIENT') 
CORRELATION  COEFFICIENT  FOR  SIGNIFICANT  • 


MEAN 
3) 


VECTOR 


VECTOR  TWO' ) 


SUBROUTINE  REL! PHI  1, PHI  2, PHI  ) 
PHI=PHI1-PHI2 

IF!  PHULT.O.C)  PHI=PHH-360. 
RETURN 
END 


SUBROUTINE  OISC  ( XBAP , M , D , PS , C ) 

D I  ME  MS  I  ON  X  B  AR  (  1  )  ,  P  S  (  1  )  ,  C  (  1  )  ,  0  (  1  ) 

SUM=0. 

DO  76  J=l, M 

MC=J+M 

76  PS!J)  =XBAR( J)-XBAR(MC) 
IJ  =  C 

DO  7  7  1=1, M 
C( I  )=0. 
DO  77  1=1, M 
IJ=IJ+1 

77  C(I  )=C (  I)  +  PS(L)*D(  IJ) 
DO  7  8  K=l, M 

MC=K+M 

78  SUM=SUM+(XBAP (K)+XBAP( MC) )*C<K) 
C(M+1)=   (  SUM,/ 2.  ) 

RETURN 
END 


10 
15 


20 


SUBRCUT INE  M INV ! A , N , D, L , M ) 

DIMENSION  A(l)  ,L(1) ,M( 1  ) 

D=1.Q 

NK  =  -N 

DO  80  K=l, N 

NK=NK+N 

L(K)=K 

M!K)=K 

KK=NK+K 

BIGA=A( KK) 

DO  2C  J=K.M 

IZ=N*1 J-l) 

DO  20  I=K,N 

I J=IZ+I 

IF(  ABSCBIGAI-  ABS( A! IJ) ) ) 

BIGA  =  A(  IJ) 

L(K)=I 

M(K)=J 

CONTINUE 


15,20,20 


38 


J=L(K) 

IF(J-K)  35,35,25 
25  KI=K-N 

DO  30  1=1, N 

KI=KI+N 

HCLD=-A(KI ) 

JI=KI-K+J 

A(KI  )=A( JI  ) 
30  A(JI  )  =  HCLD 
35  I=M(K) 

IFU-K)  45,45,38 
38  JP=N*( 1-1) 

DO  40  J  =  1,N 

JK=NK+J 

J1=JP+J 

HOLD=-A(JK) 

A( JK»=A( JI ) 
40  A( J  I  )  =HOl 0 

45  IF (BIGA I  ^8,46, A8 

46  D  =  C.C 
WRITE(6,4C0) 

400  FORMAT (4CX,« SINGULAR  INPUT  TO  MINVM 

RETURN 
48  DO  5  5  1  =  1,  N 

IF(I-K)  50,55,50 
50  IK=NK+I 

A(  IK)  =  A<  IK  )/(-BIGA) 
55  CONTINUE 

DO    6  5    1=1, N 

IK=NK+I 

HOLD=A( IK) 

IJ=I-N 

DO    6  5    J=1,N 

1J=IJ+N 

If(I-K)  6C,65,60 
60  IF(J-K)  62,65,62 
62    KJ=1J-I+K 

A(  IJ  )=HGLD*A(KJ )+A(  IJ ) 
65    CONTINUE 

KJ=K-N 

DO    7  5    J=1,N 

KJ=KJ+N 

IF(J-K)  7C,75,7  0 
70  A(KJ )  =  A(KJ  J/H1GA 
75    CONTINUE 

D=D*PIGA 

A(KK)= 1*0/ BIGA 
80    CONTINUE 

K  =  N 
100    K=(K-1 ) 

IF(K)     15C,150,105 
105    I =  I  ( K  I 

IF(I-K)     120,120,108 
108    JQ=N*(K-1) 

Jk=N*(  1-1) 

DO    110    J=1,N 

JK=JC+J 

HCLD=A( JK) 

JI=JR+J 

A(JK)=-A(JI) 
110    A( JI )     =HCLD 
120    J=M(K) 

IF(J-K)     ICC, IOC, 125 
125    KI=K-N 

DC     130    I=1,N 

KI=KI+N 

HOLD=A( KI  ) 

JI  =  KI--K  +  J 

A(KI)=-A(J  I  > 
130    A( JI )     =FCLD 

GO    TO    ICO 
150    RETURN 


39 


END 


100 


120 
130 


140 


150 
160 


170 
130 


su 

DI 
MM 
DC 
DC 
N4 
L  = 
LM 
DC 
Nl 
FN 
DO 
LM 
XR 
DU 
L  = 
XB 
XB 
LM 
DO 
LL 
DO 
LL 
N2 
CM 
LL 
DC 
DO 
LL 
D( 
N4 
LI 
DC 
LL 
FN 
DC 
D( 
RE 
EN 


PR 
ME 
=  M 
1 
I  ) 
=  0 
0 
=  0 

1 

=N 

=  N 

1 

=  L 

AR 

1 

L  + 

AR 

AR 

EA 
i 

=N 

1 

=  1. 

=  1 
EA 
=  0 
1 
1 
=  L 
LL 
=  N 

1 

=  L 

=  L 
i 

I  ) 

TU 
D 


CUT  INF  DMATX  ( K , M , N , X , XBAR , D, C HE  AN ) 
NSICN  N(1),X(1),XBAR(1),D(1),CMEAN(1) 
*M 

00  1=1, MM 
=  0«C 


60  NG=1,K 

1NG) 

1 

30  J  =  l  ,M 

M+l 

(LM)=0.0 

20  1=1, Nl 

1 

(LM)=XBAR(LM)+X(L) 

(LM)=XOAP( LM) /FN 

N=LM-M 

50  1=1 ,N1 

4+ I -Ml 

AC  J=1,M 

L+Nl 

MFAN+J 

M J)=X (LL)-XBAR(N2  ) 

50  J  =  l  ,M 

50  JJ=1,M 

L+l 

l=D(LL ) +CMEAN( J  )  *'CMEAN(  J  J  ) 

A-f  N1*M 

K 

70  1  =  1,  K 

L+N( I ) 

L 

80  1=1, MM 

=  D(  I  )/PN 

RN 


C 
C 

c 
c 

C 

c 

c 


c 
c 

C 

c 


SUBROUTINE  DAT  A ( M , N ,NS , X , XD ) 

DIMENSION  N(  2)  ,X(  5  5,23,2)  ,  SOf)(  8,66)  ,SLD(  16)  ,XD(  2116)  , 
/C(2) 

READ  (5,200)  SCD 
SOD  IS.  AN  If  PUT  MATRIX  CONTAINING  DATA  DELATING 
TO  MANNING  STRENGTH  IN  THE  ST  RATING  FOR  EACH 
DESTROYER  PARTICIPATING  IN  THE  EXERCISES  . 

THE  INPUT  MATRIX  SLD  CONTAINS  SONIC  LAYER  DEPTH 
DATA  POP  EACH  EXERCISE. 
READ  (5,100)  SLD 

C( 1)=FL0AT(N(2)  )  /FLOAT(NS) 
C(2)=-FL0AT(N( 1 ) )/FLOAT(NSJ 


CONVERT  180  DEGREES  TO  RADIAN  MEASURE 

DEN=  180. /3. 141593 

L  =  0 

DO  68  K=l,2 

NK=N(K) 

DO  6  1  1=1, NK 

L=L+1 

N2=L-NS 

RE  AD (5, 102)  (X(  I , J,K) , J  =  l,  M) 

X(I,23,K)=C(K) 


40 


c 

c 
c 


CONVERT  GEOGRAPHIC  DATA  TO  RELATIVE  DATA. 

CALL  REM  X( I , 12, K)  ,X(  I,5,K) ,RB  ) 
X(  I  ,  12, K)=P3/DEN 

CALL  PFL(X(I,13,K),X(I,5,K) , OUT ) 
OUT=OUT/DEN 

X(  I,5*K  )  =  (  X( I ,6,K)**2.+X( I, l4,K)v*2.-2.*X< I ,6,K)* 
/X(  I  ,14,K)*C0S(  OUT)  )**0.5 
DO  63  KK=1,16 
Y=FLGAT (KK ) 

IF ( X{ I  ,  1,K).EQ.Y)  GO  TO  62 
GC  TO  63 

62  X( I ,22  ,K)=SLD(KK) 

63  CONTINUE 

DO  64  KD=1 ,150 

IF( (X{ I  ,  1  ,K). EQ.SODt 1,KD) ) .AND.(X(  1 , 3 , K  )  . EQ. SOD( 2 , KD ) 
/)  )  GO  TO  7 

64  CON!  IN'UE 
7  P=6. 

X(  I,  13,K)  =  0. 
DO  74  KX=3 ,3 

X(I,13,K)=X< 1,13, K)+P* SOD (KX,KD1 
74  P=P-1. 

DO  61  JJ=1 ,M 
N2=N2+."!S 
61  XD(N2)=X(I,JJ,K> 
68  L=NK 
100  FORMAT ( 7F1C.0) 
102  FORMAT  (  3F1C.  0,  1  OX,  3F1  CO,  /,7F  10,0,/,  7E  10,0,/,  F  10.0, 

/1CX,2F  ICC  ) 
200  FORMAT (F 4. C,F6.0,6F10.0) 
RETURN 
END 


100 


102 


105 


107 
108 


110 


115 

127 

130 

135 

137 


SUBROUTINE  CORPE  (N,M, 
DIMENSION  X(l) ,XBAR(1 ) , 
/T(  1) 
DO  100  J=1,M 
B( J) =0.0 
T( J)=C.O 
K= (M*M+M )/2 
DO  102  1=]  ,K 
R(  I  1^0.0 
FN=N 
L  =  0 

IF(IO)  105,  127,  105 
DO  ice  J=1,M 
DO  107  1=1, N 
L  =  L  +  1 

T(J)=T(J)+X(L) 
XBARi J) =  T< J) 
T( J)=T( J  )/FN 
DO    lit     1=1 ,N 
JK  =  0 
L=I-N 

DO    11C     J=l »M 
L  =  L  +  N 

D( J)=X(L)-T(J) 
B( J)=B( J)+ 
DO  115  J=l 
DO  115  K=l 
JK=JK+1 

R( JK  >  =  R  ( JK  )+D( J  )-D(K) 
GO  TO  2  05 

IF(N-M)  13 C,  130,  135 
KK  =  N 

GC  TO  137 
KK=M 

DO  14C  1=1 ,KK 
CALL  DATA  (M,D) 


10, X, XBAR,STD,RX,P 
STD( 1 ),PX(1  ),R(  1  ) 


,B,D,T) 

,  B  (  1  )  ,  D  (  1  )  , 


C)<  J) 

»M 

,J 
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140 
150 

170 

180 
185 


190 


200 
205 


210 


220 


DC  14C  J  =  l  ,M 

T( J)=T( J)+D( J) 

L=L  +  1 

P  X  (  L  )  =  J  <  J ) 

FKK=KK 

DO  150  J=l ,  M 

XBAP( J) =T( J) 

T{ J  >  =  T( J)/FKK 

L=0 

DO  160  1=1 iKK 

JK=0 

DO  17C  J=l ,M 

L  =  L+1 

D(J)=PX(L)-T(J 

DO  180  J  =  l »M 

B( J)=B( J )+D( J) 

DO  18C  K  =  l  ,J 

JK=JK+1 

R( JK  )= >  ( JK  )+D( 

IF(N-KK)  2  05, 

KK=N-KK 

DO  200  1=1 , KK 

JK  =  0 

CALL  DATA 

DC  190  J=l 


222 

225 

230 


240 


2  50 


) 


J  )«D(K) 

205,  185 


ABS(RUK)  )  ) 


(M,D) 

»M 
XBAK( J) =XPAR( J )+D( J) 
D(J)=D(J)-T(JJ 
B(J)=B(J)+D(J) 
DO  2  00  J  =  l  ,M 
DO  200  K  =  l  ,  J 
JK=JK+1 

R{ JK)=R(JK)+D( J)*D(K) 
JK  =  C 

DO    2  1C     J=1,M 
XBAf<  (J)  =XBAR( J)  /FN 
DO    210    K  =  l  ,J 
JK=JK+1 

F  (  JK  )  =  R  <  JK  )  -  B  (  J  )  *  B  (  K  )  /  F  N 
JK  =  C 

DO  220  J=l ,M 
JK=JK+J 
STD(J)=  SQRT( 
DO  2  30  J=l  fM 
DO  230  K  =  J  ,M 
JK=J+( K*  K-K  )/2 
L=M*(J-1)+K 
RX(L  )  =  P ( JK  ) 
L=M*(K-1 )+J 
RX(LJ=R( JK) 

IF(STD( J)*  -STD(K)  )  225,  222, 
R (JK )=0.C 
GO  TO  2  30 

R(JK)=P(JK)/{STD(J)*STO(K)) 
CONTINUE 
FN=SQRT(FN-1.0) 
DC  2*C  J  =  l  ,M 
STD( J)=STDU)/FN 
L  =  -M 

DO  250  1=1, M 
L=L+M+1 
B  (  I  )  =R  X  (  L ) 
RETURN 
END 


225 


SUBROUTINE     ORDER    {  M  ,  R  ,  NDP  P,  K  ,  I  S  AVE  ,  RX  ,RY  ) 
DIMFNSION    P( 1)  ,  ISAVE( 1 ) ,RX( 1  ),RY(1  ) 
MM  =  0 

DO    13C     J=1,K 
L2  =  ISAVE(J  ) 

IF(NDEP-L2)     12  2,     123,     123 
122    L=NDCP+<L2*L2-L2)/2 
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123 
125 


127 

128 
129 
130 


GG  TO  125 

L=L2+(NDEP*NDEP-NDEP)/2 

RY  ( J  )  =  R (  L) 

DO  13C  1=1  fK 

L1=ISAVE(I ) 

IF(L1-L2)  127,  123,  128 

L=L1+(L2*L2-L2 

GO  TO  129 

L=L2+(L1*L1-L1 

MM=MM+1 

PX(MM»=R(L ) 

ISAVE(K-H)  =NDEP 

RETURN 

END 


123, 
)/2 

)/2 


,STD(1 
SR(1  ) 


SUBROUTINE 
/ANS) 

DIMENSION    XBAR(l) 
/T( 1)  ,AMS(1  ) ,3<  1)  , 

MM=K+1 

DO    100    J  =  l  ,K 
100    B(J)=0.f 

DC    110    J  =  l  ,K 

L1=K*{ J-l) 

DO    11C     1=1  ,K 

L=L1+I 
110    B( J)=E( J)*KY( I )*RX(L> 

RM=0.0 

BO=C.O 

L  1=1  SAVE  (f> 

DO    120     1=1 

RM=RM+B< I )^RY(  I  ) 

L=ISAVE (  I) 

IH1)  =  B(I)*(STD(L1)/STD(LI) 
120    BO=BO+B<  I  )  >  XbAR(L) 

BG=XBAR(L1  )-B0 

SSAF=RM*D(L1  ) 
122    RM=    SOf-T(     ABS(RM)  I 

SSDR=D(L1)-SSAR 

FN=M-K-1 

SY=SSDR/FN 

DO     13C    J  =  l  ,K 

L1=K*(J-1)+J 

L=ISAVE( J) 

Sb(  J  )=     SOFT(     ABS( 

T( J>=B( J)/S3( J  ) 


MULTR     (N,K,XBAR,STD,D,RX,RY,ISAVc,B,Sb,T 


) ,D(1  ),RX(1  ),RY(1  )  ,  ISAVE(  1  ) 


125 
130 
135 


M) 


(RX(L1)/D( L) )*SY) ) 


SY=     SQRTI     A3S(SY)) 

FK=K 

SSAPM=SSAF/FK 

SSDRM=SSDP/FN 

F=SSAPK/SSDRM 

ANS(1)=B0 

ANS( 2)=RM 

ANS( ?)=SY 

ANS(*)=SSAF 

ANS(5)=FK 

ANS(6)=SSAPM 

ANS(7)=SSDR 

ANS( 8)=FN 

ANS(9)=SSDRM 

ANS( 1C^=F 

RETURN 

END 
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